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Abstract 

We continue the analysis of S-wave production amplitudes for the 
reaction 7r~p — > 7r + 7r~n involving the data obtained by the CERN- 
Cracow-Munich collaboration on a transversely polarized target at 17.2 
GeV/c 7r~ momentum. This study deals with the region below the KK 
threshold. In particular, we study the "up-steep" solution containing 
a narrow S-wave resonance under the p(770) . This solution exhibits a 
considerable inelasticity rj which does not have any physical interpreta- 
tion. Assuming that this inelasticity behaviour represents an unlikely 
fluctuation we impose rj = 1 for all data points. This leads to non- 
physical results in one third of the 7r + 7r _ effective mass bins and in the 
remaining mass bins some parameters behave in a queer way. The sit- 
uation is even worse for the "down-steep" solution. We conclude that 
the 17.2 GeV data cannot be described by a relatively narrow /q(750). 
The "down-flat" and "up-flat" solutions which easily pass the 77 = 1 
constraint exhibit a slow increase of phase shifts in the p(770) mass 
range. 

1 Introduction 

Scalar mesons are one of the main puzzles of light quark spectroscopy. 
Even in the lowest mass region (below the KK threshold) the situation is far 
from being clear. In addition to a broad /q(500) interpreted either as a q — q 
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object (see e.g. Ref. |T|) or as a glueball by Ochs @, the relatively narrow 
/o(750) has been persistently claimed by Svec [3-5] . Arguments against the 
narrow / (750) has been given e.g. by Morgan in Ref. ||. The main source of 
information in this mass region is the irn partial wave analysis (PWA) yielding 
the S"-wave. It should be stressed that a study of >S-wave objects does require 
the partial wave analysis to "subtract" the dominant contribution of leading 
p(770) meson. 

Virtually all PWA's in last decades were based on the old CERN-Munich 
experiment 0], which supplied 3xl0 5 events of the reaction 

7T~p — > 7T + 7T~n (1) 

at 17.2 GeV/c. The number of observables provided by such experiment is 
much smaller than the number of real parameters needed to describe the partial 
waves. Consequently, the dominance of pseudoscalar exchange, equivalent to 
the absence of pseudovector exchange and several other physical assumptions 
have been made in previous studies [7-8] . These results have been generally 
used without even mentioning the assumptions essential for their derivation. 

In our previous paper || (hereafter called paper I) we have used the results 
of PWA performed in the the effective mass m n7T range from 600 MeV to 1600 
MeV at four-momentum transfer squared \t\ = (0.005 — 0.200) GeV 2 /c 2 using 
additionally the results of the polarized target experiment. This experiment, 
performed 25 years ago by the CCM (CERN-Cracow-Munich) collaboration 
[110], provided 1.2xl0 6 events of the reaction 

7T~Pf — > 7T + 7T~72 (2) 

also at 17.2 GeV/c. Combination of results of both experiments yields a num- 
ber of observables sufficient for performing a quasi-complete and energy inde- 
pendent PWA. This analysis is only quasi-complete because of an unknown 
phase between two sets of transversity amplitudes. Nevertheless, full (contain- 
ing both 7r and a x exchange) intensities of partial waves could be determined in 



a model-independent way. The original study of the CCM collaboration |TT 
removed ambiguities appearing in earlier studies, except for the "up-down" 
ambiguity JTSJ. The "up" solution contains an 5*-wave resonance just under 
the p(770) and of similar width, while the "down" S*-wave modulus stays high 
and nearly constant all the way to the /o(980) . 

In paper I we have made a further step in the analysis of 17.2 GeV/c 
data bridging two sets of transversity amplitudes. We required the phases 
of the leading P, D and F-transversity amplitudes to follow the phases of 
the Breit-Wigner p(770) , ^(1270) and ps(1690) resonant amplitudes in the 
low, medium and high mass region, respectively. Further, using the measured 
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phase differences between the S-wave and the higher waves we determined 
the absolute phases of the S-wave transversity amplitudes. Once the phases 
are known, the S'-wave amplitudes of different transversity can be combined 
which allows us to determine explicitly for the first time the pseudoscalar and 
pseudovector exchange amplitudes in the S'-wave. This has been done using 
much weaker assumptions^ than those made in any earlier analysis. The price 
we pay is a fourfold ambiguity in our pseudoscalar exchange S-wave ampli- 
tude. In addition to "up-down" ambiguity of the old CCM analysis |l]J there 



are ambiguities resulting from adding or subtracting the phase difference since 
the PWA yields only the absolute value of the phase difference. Thus we have 
"down-flat", "down-steep", "up-flat" and "up-steep" solutions. Differences be- 
tween "flat" and "steep" refer mainly to the behaviour of the S-wave phase 
shifts below the /o(980) . Above the /o(980) all the solutions are fairly simi- 
lar. It is the region below the / (980) which is a subject of the present paper. 
The main difference is that both "steep" solutions contain a relatively narrow^] 
resonance under the p(770) (like the old "up" solution) while both "flat" solu- 
tions indicate a /o(500) state with a width of about 500 MeV. In particular the 
"down-flat" solution is very similar to the old solution of the CERN-Munich 
group 0. 

In paper I we have determined inelasticities rj for isoscalar 7r-exchange am- 
plitudes in all solutions (see Fig. 1). It is obvious that both "flat" solutions 
easily pass the r\ < 1 test and the "down-steep" solution is not acceptable. 
Unfortunately, the situation is not as simple as presented in the last edition 
of Review of Particle Properties [|T^]. The authors of "Note on scalar mesons" 
discussing a hypothetical narrow state at 750 MeV write: " Such a solution is 
also found by (KAMINSKI 97)... However they show that unitarity is violated 
for this solution; therefore a narrow light fo state below 900 MeV seems to be 
excluded". The point is that our "up-steep" solution although exhibiting "a 
puzzling behaviour of inelasticity" cannot be excluded so simply as our "down- 
steep" solution and it also contains a narrow / (750). It should be recalled, 
that contrary to the Svec analysis which uses only moduli of the unseparated 
(pseudoscalar and pseudovector exchange) S-wave, we study inelasticity and 
phase shift of the pure itn — > tttt isoscalar S-wave. In this paper we discuss the 
feasibility of an interpretation of 17.2 GeV data in terms of a narrow /o(750). 

The paper is organized as follows. In Sect. 2 we study the inelasticity 
behaviour in more detail. In Sect. 3 we impose strict r\ = 1 condition on all 
solutions for the effective nn masses below the KK threshold. The results are 
discussed in Sect. 4 and summarized in Sect. 5. 



Hhe main assumption was neglecting of a possible influence of the a\ exchange in / 
S'-wave as well as in P, D and F waves. 

2 hereafter "relatively narrow" means "with a width close to T p = 150 MeV". 



3 



2 5- wave inelasticity 



As seen in Fig. la the inelasticity for the "up-steep" solution behaves in 
a non-trivial way. Below 720 MeV and above 820 MeV rj > 1 while in the 
intermediate mass region rj < 1. Qualitatively this behaviour of the "up-steep" 
solution is very similar to the "down-steep" solution. The latter solution was 
excluded in paper I since the corresponding inelasticity significantly exceeded 
unity for m n7V > 820 MeV. A general trend of 77 points for the " up-steep" solu- 
tion resembles very much a trend of the "down-steep" data in the region below 
820 MeV. Above 820 MeV the "up-steep" values of inelasticity lie also above 
unity. They are, however, closer to unity than the "down-steep" inelasticities. 
This fact prevented us to reject in paper I the "up-steep" solution solely on a 
basis of the minimization fit of the sum ~ I) 2 i n the whole energy region 
between 600 MeV and 1000 MeV. 

Now we examine in detail the m nn dependence of the inelasticity i] corre- 
sponding to the four solutions. Let us define by a, b and c three ranges of 
WW 600 < < 720 MeV, 720 < m n7T < 820 MeV and 820 < m W7V < 1000 
MeV, respectively. Then we make a simple fit to inelasticities in the above 
three ranges by a constant N. We minimize the sum J2i(Vi ~ N) 2 /Ar]i 2 , where 
Ar/i are the experimental errors of rj. The results are shown in Table 1. Let 
us notice that for the "steep" solutions the second value in the intermediate 
region b is definitely lower than unity, this is also different from the values in 
the ranges a and c which in turn are higher than unity. In this aspect the 
"up-steep" solution is very similar to the solution "down-steep" already re- 
jected in paper I. On the other hand both solutions "down-flat" and "up-flat" 
do not exhibit any dip of 7] in the range b; the constants N are very close to 
unity everywhere. If we determine one common value of N in the whole range 
between 600 MeV and 1000 MeV then we obtain 1.00±0.06 for the solution 
"down-flat" and 0.98±0.06 for the solution "up-flat". Both are compatible 
with unity as seen in Fig. lb. 

Table 1: Constants N fitted to rj in three ranges: 600 < m ww < 720 MeV (a), 
720 < m W7V < 820 MeV (b) and 820 < m ww < 1000 MeV (c) 



Solution 


a 


b 


c 


" up-steep" 
" down-steep" 
"up-flat" 
"down-flat" 


1.17 ± 0.10 
1.40 ± 0.12 
0.97 ± 0.10 
0.97 ± 0.13 


0.67 ± 0.17 
0.82 ± 0.10 
1.02 ± 0.11 
1.05 ± 0.09 


1.18 ± 0.10 
1.67 ± 0.11 
0.95 ± 0.12 
0.96 ± 0.09 
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Figure 1: a) Scalar-isoscalar irn inelasticity coefficient i] versus the effective nir 
mass for the "down-steep" (full circles) and "up-steep" (open circles, shifted 
by 6 MeV) solutions. Solid and dashed lines represent fits of the fourth order 
polynomial to the coefficient 77 for the "down-steep" and "up-steep" solutions, 
respectively, b) Same as in a) but for the "down-flat" (full circles) and "up- 
flat" (open circles) solutions. Solid line shows the rj — 1 value fitting the 
experimental data. 5 



Obviously we cannot obtain constant fits to rj of a similar quality for the 
"steep" solutions as those fits to the "flat" solutions. In order to probe rather 
strong energy variation of 77 in the former solutions we have tried to fit it 
by polynomials of different powers. We treat these fits only as an ad hoc 
description of data in the particular region between 600 and 1000 MeV. It 
turned out that in order to obtain good fits we need polynomials of the order 
as high as four (see Fig. la). Once again one can see the similar dependence 
of inelasticities when one compares a shape of the fitted curves to the "up- 
steep" and "down-steep" data points. A minimum of r\ for m n7r around 800 
MeV, indicating an inelasticity, will be discussed below. Before passing to a 
physical discussion of possible consequences of the strong energy dependence of 
inelasticity for the "up-steep" solution below the KK threshold, let us discuss 
a possibility of the fluctuation of rj values around unity. It is true that the 
deviation is not significant (see paper I: x 2 — 15 for 17 points below 940 MeV) 
if we ignore the shape of the 77 distribution. However, the chance of all five 
7] < 1 points falling accidentally in the effective mass region between 720 and 
820 MeV is only 2 ■ 10~ 3 . This value was calculated by taking into account a 
number of all the possibilities to choose five lowest values of 77 among twenty 
available points in the effective mass range between 600 and 1000 MeV and 
then grouping them together in the range b. 

In paper I in addition to inelasticities we have calculated the phase shift 
values. For both "flat" solutions phase shifts grow slowly with m^. For the 
"steep" solutions we have obtained rather fast increase of the S-wave 1=0 tttt 
phase shifts near 770 MeV - the value close to the p-meson mass. We have 
tentatively fitted inelasticities and phase shifts of the solution "up-steep" by a 
single resonance. This was done for m n7T < 940 MeV in order to avoid a possi- 
ble influence of the /o(980) . We additionally allowed a simultaneous change of 
all 77 values by the same factor R v since in paper I they were fixed only by min- 
imizing the Sj(?7j — l) 2 value. In fact a combined fit (x 2 /NDF = 26/30) yields 
the reduction factor R n = 0.71 ± 0.10 and the following resonance parameters: 
m = (754 ± 5) MeV, T = (162 ± 9) MeV and x = 0.74±g;Jg, x being an elas- 
ticity of the resonance. The mass and the width agree with m = (753 ± 19) 
MeV, T = (108 ± 53) MeV claimed by Svec ||, on the basis of the same data. 
However, such a considerable inelasticity is inconsistent with the available ex- 
perimental data on the 4tt system which is the only kinematically possible 
channel. 

The lowest mass of the 47r system is probably available in the central pro- 
duction due to a l/m 2 ^ flux factor. In fact the WA91 [111 and WA102 ITl 



collaborations have found a tiny peak around 800 MeV in the mass distribu- 
tion of their Arc system produced at 450 GeV in the reaction 

PP — > PfTY + 'K + 'K~'K~p s , (3) 
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where pf and p s stand for fast and slow proton, respectively. This is however 
well explained by the reflection from the 77 — > i]7i + 7i~ decay with the loss of 
the slow 7T° from the subsequent 77 — > 7r + 7r~7r° decay. The IHEP-IISN-LANL- 
LAPP collaboration [|T^] studying at low \t\ a reaction similar to reaction (1) 
i.e. 

n~p -> TrVVVVi (4) 

has found that the effective mass distribution of their An system starts only 
around 800 MeV and rises smoothly. The LRL group (T7j studying the reaction 



n + p -> 7T + 7r + 7r _ 7r _ A ++ (5) 

found hardly any events with a 4n mass below 1 GeV. The 4-7T mass spectrum 
from annihilation pN — > 5n starts at even higher mass as shown in [18-21] . 
Thus the non-zero inelasticity in the "up-steep" solution does not have any 
reasonable physical interpretation. In next section we will assume 77 = 1 . 



Another approach to the tttt scalar-isoscalar 
phase shifts 



In Sect. 2 we have shown that the Air channel is very weak below 1000 
MeV. Therefore now we assume that the titt S-wave inelasticity is exactly 
equal to unity up to 990 MeV and we shall make a new analysis of the Tin 
isoscalar-scalar phase shifts obtained from the ir~p — > -K + i{~n data at 17.2 
GeV/c. Let us recall that in paper I a separation of the S-wave pseudoscalar 
A and pseudovector B exchange amplitudes for the production process was 
performed and that we have calculated the S-wave tt + 7t~ — > tt + 7t~ amplitude 
as from the following formula: 

a s = KfA . (6) 
In (§) K is the proportionality factor 

%-npTty/sq; 1 

where p n is the incoming 7r~ momentum in the ir~p centre of mass frame, s 
is the square of the total energy in the same frame, q n is the final pion mo- 
mentum in the irir rest frame, |- is the pion-nucleon coupling constant (taken 
as 14.6 in paper I) and M is the proton mass. The coefficient / is the com- 
plex correction factor which (when averaged over the four momentum transfer 
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squared) represents the t-dependence of the pion-nucleon vertex function, the 
off-shellness of the exchanged pion and a possible phase of the pion-exchange 
propagator — a case where at high energy the exchanged pion is treated as a 
Regge particle. This coefficient was calculated from the requirement that the 
sum J2i{Vi ~ I) 2 was minimal for the inelasticities of the scalar-isoscalar txti 
amplitude ao for a set of points depending on the tttt effective mass up to the 
KK threshold mass. The amplitude a is connected to the amplitude a s and 
the isospin 2 S-wave amplitude a 2 in the following way: 

1 

a = 3a s - ~a 2 . (8) 
The ao is also related to the isospin S-wave phase shift Sq and inelasticity 77: 

ne 2t5 ° - 1 

* = . 0) 

Motivated by the results of Sect. 2 we impose now the condition rj = 1 in the 
whole range below 1 GeV. Then, from (S) 

a = sin5oe i(50 , (10) 

so the modulus of a is uniquely related to the phase shift value 5q. Since 
in (HP A is the complex amplitude calculated with some errors coming from 
experimental uncertainties, then (^) and fllPP are not necessarily satisfied if we 
require rj = 1. In order to keep this assumption valid we have to include at 
least one additional real factor n in equation (g|) at each m nn , so now 

a n s ew = na s . (11) 

We shall also assume that the isospin 2 amplitude a 2 is fully elastic and 
can be described by the corresponding phase shift 82 

a 2 = sm8 2 e iS2 . (12) 
Then, following (0) and inserting (JET]) into (^) one has to satisfy equations 

1 2 



2 1 1 2 

77 = |1 + 2zao\ 



1 + 2t(3KfnA - ~a 2 ] 



1. (13) 



We treat fll3[) as a quadratic equation for n. Its roots are 

1 

n = — 
6 

where 



n = —^-(b±Jb 2 -3sm 2 S 2 ), (14) 
6|a s | 

b = (1 + sin 2 5 2 ) sin a + - sin 2 <5 2 cos a (15) 
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and a denotes the phase of as 



a s =\a s \e ia . (16) 

Let us remark that in the limit of 82 going to (vanishing 02) we should 
consider only the upper sign (+) in fll^l) . In general, both solutions for n are 
possible; we have, however, used only the root with the upper sign since its 
value was closer to unity. The isotensor phase shift 82 is calculated according 
to the parametrization given in paper I which fits well the data of Ref. |22| 
obtained by method B. 

For each value of m nn we have to check whether the roots exist. With rj = 1 
the amplitude a must satisfy the elastic unitarity condition 

Im a = \a \ 2 . (17) 

Therefore the following inequality must be fulfilled 

\b\ > V3|sin5 2 |. (18) 

We have obtained the following numerical results for the solutions discussed 
in I: the inequality (|18|) was satisfied for all twenty points of the solution 
"up-flat" and for 19 points (except of the extreme point at 990 MeV) corre- 
sponding to the solution "down-flat" . However for 7 points of the solution 
"up-steep" and for 12 points of the solution "down-steep" the condition ( |1"8"D 
was violated. This fact casts a serious doubt on a validity of both "steep" 
solutions. The resulting values of n calculated in cases when ( |T8| ) was satisfied 
are shown in Fig. 2. The errors of n are due to experimental errors of the 
modulus I a,s I and the phase a extracted from experiment. No errors of 82 were 
taken into account since 82 values were calculated using the smooth theoretical 
parametrization. We see in Fig. 2b that both "flat" solutions are well fitted by 
constants very close to unity (0.994±0.03 for the "down-flat" and 0.997±0.04 
for the "up-flat" solution). On the other hand a variation of n with m n7T for 
two "steep" solutions is better described by a parabola than by a constant (see 
Fig. 2a). Such strong dependence on of the coefficient n corresponding 
to both "steep" solutions can be used as a fairly strong argument against an 
acceptance of these solutions as good physical solutions. 

For completeness we present in Fig. 3 new irn phase shifts calculated from 



( |10"D using alf w given by ([TTJ) . Obviously we show only these points for which 
the corresponding values of n do exist. The new phase shifts, calculated under 
the assumption that rj = 1, agree very well with those presented in paper I for 
the "flat" solutions. For the "steep" solutions this agreement is not so good 
and the new errors of 8q are larger than those shown in paper I. 
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Figure 2: a) Effective mass dependence of the parameter n for the "down- 
steep" (full circles) and "up-steep" (open circles, shifted by 6 MeV) solutions. 
Solid and dashed lines represent fits of the second order polynomial to n for 
the "down-steep" and "up-steep" solutions, respectively. Note that in many 
bins no physical solution could be found, b) Same as in a) but for the "down- 
flat" (full circles) and "up-flat" (open circles) solutions. Solid line represents 
values of the constant parameters fitted to n for the " down-flat" and " up-flat" 
solutions. i n 
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Figure 3: Comparison of the scalar-isoscalar nn phase shifts obtained in pa- 
per I (circles, shifted in m n7T by 6 MeV) with the new phase shifts (squares) 
calculated under the assumption r] = 1. a) For the "up-steep" solution. Solid 
line shows the Breit-Wigner fit to the data marked by open circles as described 
in Sect. 2. b) For the "down-steep" solution, c) For the "up-flat" solution, 
d) For the "down-flat" solution. 
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4 Discussion 



In Sections 2 and 3 we have presented arguments that both "steep" so- 
lutions have unphysical behaviour. On the other hand both "flat" solutions 
satisfy well our tests and none of them can be eliminated using the methods 
described in this paper. Therefore let us discuss common features of these 
solutions and major differences between them. In Fig. 4 the "flat" solutions 
are plotted in a wide effective mass range up to 1600 MeV. Their shape is 
quite similar. One sees an initial steady grow of phase shifts with m n7T above 
600 MeV, then at about 1000 MeV, corresponding to the KK threshold, there 
is a jump as high as 140° and further on a fairly steep increase above 1300 
MeV. An interpretation of this behaviour of phases in terms of three scalar 
resonances /o(500), /o(980) and / (1500) , started in paper I, has been contin- 
ued in more detail in Refs. [23,24] . We do not repeat it here but we underline 
major differences between the "up-flat" and "down-flat" phase shifts since they 
lead to different values of the /o(500) resonance parameters. The /o(500) mass 
for the "up-flat" solution is by about 50 MeV higher than the corresponding 
mass for the "down-flat" solution. The reversed relation for the /o(500) width 
leading to a difference between 45 and 50 MeV is also observed. We do not 
see important differences between the /o(980) parameters for the above solu- 
tions. As seen in Fig. 4 the most important differences, reaching about 45°, 
exist between 800 and 1000 MeV. They are related to a difference between the 
moduli of the S-wave pseudoscalar amplitude (see paper I). This difference of 
moduli is then directly transformed into a difference between the phase shifts 
because inelasticity for both "flat" solutions below the KK threshold is very 
close to unity. This fact in turn guarantees a fulfilment of the elastic unitarity 
condition (|i~7l) . Closer inspection into Fig. 4 allows one to see a continuity 
of the phase differences above the KK threshold, namely the " up-flat" points 
lie systematically above the "down-flat" points up to about 1300 MeV. Above 
this value of we do not observe any systematic difference. 

Since the data points of the "up-flat" solution lie between the points of 
the (already excluded) "up-steep" solution and the "down-flat" one, we have 
checked whether the /o(750) survives in this solution. This was done by fitting 
elasticities and phase shifts of the "up-flat" solution by a single Breit-Wigner 
term like it was done in Sect. 2 for the "up-steep" solution. No overall change 
of r] values was needed (R v = l.Oll&iI) 

and the resonance parameters are 
m = (732 ± 8) MeV, T = (2461H) MeV and x = l.OOlg;^. The large width is 
inconsistent with a narrow / (750). 
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Figure 4: Comparison of the scalar-isoscalar nn phase shifts obtained in 
paper I for the "down-flat" (full circles) and "up-flat" (open circles) solutions. 
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5 Summary 



In conclusion, we have studied in detail the tttt effective mass dependence 
of the S-wave isoscalar phase shifts corresponding to four solutions "up-steep", 
"down-steep", "up-flat" and "down- flat" found in paper I. Both "steep" so- 
lutions exhibit an inelasticity behaviour which has no physical interpretation. 
We do not find any data on the 4ir systems which could explain a strong m n7T 
dependence of the inelasticity corresponding to the "steep" solutions below 
the KK threshold. The "down-steep" solution was already rejected in paper 
I since its inelasticity substantially exceeded unity for above 820 MeV. 
Assuming that the "up-steep" inelasticity is an unusual fluctuation we impose 
7] = 1 for all points and, for completeness, in all solutions. This leads to non- 
physical results in 7/20 mass bins for the "up-steep" solution and in 12/20 
bins for the "down-steep" solution. In the remaining mass bins the parame- 
ters behave in a queer way (compare Figs. 2a and 2b ). We conclude that the 
"up-steep" solution cannot be treated as a good physical set of phase shifts. 
It can be eliminated together with the "down-steep" solution. However, the 
"up-flat" and "down-flat" solutions easily pass our tests. We would like to 
stress that both the / (500) and / (980) resonances are present in the "flat" 
solutions. This is not true for a relatively narrow / (750) . 
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